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We give a simple and direct treatment of insensitivity in stochastic networks which 
is quite general and which provides probabilistic insight into the phenomenon. In the 
. . . case of multi-class networks, the results generalise those of Bonald and Proutiere (2002, 

^ ! 2003). 
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^ ; 1 Introduction 
O , 

■ It is well-known that many stochastic networks — notably queueing and loss networks — 

■ have stationary distributions of their level of occupancy which depend on certain input 
distributions only through the means of the latter. This phenomenon of insensitivity has 

• been studied by various authors over an extended period of time, in varying degrees of 
p ^ \ generality and abstraction, and using a variety of techniques. 

r~| . In the present paper we revisit this topic to develop an insight of Pechinkin (1983, 1987) 

, to give a very simple and direct treatment of insensitivity. In particular the approach 

avoids those based on brute-force calculations, the consideration of phase-type distribu- 
tions (Schassberger, 1978, Whittle, 1985, Bonald and Proutiere, 2002, 2003), or the use of 
quite complex machinery for handling generalised semi-Markov processes (Burman, 1981, 

• Schassberger, 1986) — although such processes are implicit in the current approach. It fur- 
CN ■ ther avoids assumptions about, for example, continuity of distributions, necessary for some 

of the above approaches, and also explicitly identifies the entire stationary distributions 
of the networks concerned, showing that, where insensitivity obtains, these stationary dis- 
tributions have a particularly simple and natural form. Pechinkin used his insight, which 
involves what is in effect a coupling argument together with induction, to give probabilistic 
proofs of the insensitivity of a number of single-class loss systems with state-dependent ar- 
rival rates — results originally proved analytically by Sevastyanov (1957). He also indicated 
the wider applicability of the approach in the single-class case. In the present paper we 
give a substantial reformulation of the underlying idea, under more general conditions and 
showing that its most natural expression is in terms of balance equations. This consider- 
i ably simplifies its application to single-class systems — notably the quite complex coupling 

constructions are no longer needed. It further makes possible the extension of the idea to 
the multi-class networks considered in Section [3j The main aim is to provide probabilistic 
insight, notably for multi-class networks. Indeed it is shown that insensitivity is simply a 
byproduct, under appropriate conditions, of probabilistic independence. 

We study networks in which individuals arrive at various classes at rates which may 
depend on the state of the entire system, bringing workloads which are independent and 
identically distributed within classes and which have finite means. Within each class 
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workloads are reduced at rates which may again be state-dependent (when the rate is 
constant workloads may be identified with lifetimes in classes), and on completion of its 
workload an individual moves to a different class or leaves the system, with probabilities 
which may yet again be state-dependent. 

In order to obtain insensitivity we typically require that an individual joining a class is 
immediately served, i.e. has its workload reduced, at a rate which is the same as that 
of an individual immediately prior to leaving the class (where in each of these cases the 
number of individuals in each class of the system is the same) — more generally that the 
service discipline should define a network which is symmetric in the sense of Kelly (1979). 
The most common example is that of processor-sharing networks, but other possibilities 
are well-known, for example, "last-in- first-out preemptive resume" networks. We shall 
concentrate on a very broad class of processor-sharing networks, introduced by Bonald 
and Proutiere (2002) and including, for example, traditional loss networks and processor- 
sharing Whittle and Jackson networks, as special cases). We shall also indicate the simple 
modifications required to deal with other possibilities. 

For the above class of processor-sharing networks, Bonald and Proutiere used phase-type 
arguments to show that, under conditions which correspond to the satisfaction of the 
appropriate partial balance equations, the stationary distribution of the number of indi- 
viduals in each class is insensitive to the workload distributions, subject to the means of 
the latter being fixed and to the distributions themselves being drawn from the broad class 
of Cox-type distributions (dense in the class of all distributions on M+). In the present pa- 
per we formally consider all workload distributions on M-|_ with finite means, and identify 
also the stationary residual workload distributions. However, as stated above our main 
aim is to give a direct and probabilistically natural treatment. It turns out (and is in 
many cases well-known) that, when the appropriate partial balance equations are satis- 
fied for such a network, then the stationary distribution of the entire system, including 
the specification of residual workloads, is such that departures from each class are exactly 
balanced by arrivals to that class — in a sense again to be made precise below. Indeed, for 
single-class systems, this is the essence of Pechinkin's insight. What is of interest is that 
same idea extends to establish insensitivity for the very much more general networks con- 
sidered here, and indeed appears also to establish insensitivity in more abstract settings 
such as that considered by Whittle (1985), though we do not formally consider this more 
abstract environment here. 

In order to fix ideas, it is convenient to consider first, in Section [2l single-class networks. 
Here the extension of previous ideas is not too difficult. Nevertheless it is desirable to give 
a careful treatment of this case, avoiding notational complexity while preserving rigour, so 
as both to establish the underlying principle and also to set the scene for the multi-class 
networks which we consider in Section O 

2 Single-class networks 

Consider an open system with a single class of individual (customer, call, or job). Indi- 
viduals arrive as a Poisson process with state-dependent rate a{n), where n is the number 
of individuals currently in the system. Arriving individuals have workloads which are 
independent of each other and of the arrivals process with a common distribution fj, on 
which we assume to have a finite mean m{iJ,). While there are n individuals in the 
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system, their total workload is reduced at a rate /?(n) > 0, where we assume /?(n) > if 
and only if n > 0; an individual departs the system when its workload is reduced to zero. 
By suitably redefining the rates /3(n) if necessary, we may, and do, assume without loss of 
generality that the mean workload m(/x) = 1. 

We consider first the processor-sharing case. Here when there are n > individuals in the 
system, the workload of each is simultaneously reduced at a rate /3(n)/n, and the set-up 
described above becomes a fairly general description of a single-class processor sharing 
system. A special case is the simple Erlang loss system, in which, for some a,P>0, we 
have a(n) = al{n < C) for some capacity C < oo (where I is the indicator function) and 
/3(n) = n/3 for all n > 0. Here individuals are typically referred to as calls, and workloads 
correspond to call durations (since (3{n)/n is independent of n). A further special case 
is the M/GI /m/oo processor-sharing queue, in which, again for some a,/? > 0, we have 
a{n) = a for all n and (5(n) = min(n, m)P for all n and some fixed m. 
We represent the system as a Markov process {X{t))t>o by defining its state at any time t 
to be the number n of individuals then in the system together with their residual workloads 
at that time. (An alternative is to record, for each individual, the workload completed at 
time t.) For given n > these workloads form an (unordered) set, and may be regarded as 
taking values in the quotient space Sn obtained from M" by identifying points which may 
be obtained from each other under permutation of their coordinates. The cr-algebra B{Sn) 
on Sn is similarly formed in the obvious manner from the Borel cr-algebra on M". The 
state space S for the process {X{t))t>Q is then the union of the n > 0, where the set 
is taken to consist of a single point, and its associated o"- algebra B{S) consists of those sets 
which are countable unions of sets in the cj-algcbras B{Sn)- The process {X{t))t>Q is thus 
an instance of a piecewise-deterministic Markov process (Davis, 1984, 1993). However, we 
avoid the need for most of the general machinery for handling such processes. 
We define the probability distribution Jx on M+ to be the stationary residual life distribu- 
tion of the renewal process with inter-event distribution /i, that is, if has distribution 
function F then /i has distribution function G given by 



(recall that m{^) = 1). Note that the "residual life" here should be thought of as a residual 
workload rather than a time. For each n > 1, define also the probability distribution /i„ on 
Sn to be the product of n copies of the distribution /2, again with the above identification of 
points in W]^ (more formally, /x„(^) = /2"(^^^(A)), A G B{Sn), where /2" is the product of 
n copies of the distribution p, and 6 is the projection from into the quotient space Sn)- 
Thus fin represents the joint distribution of the residual lives at any time in a set of n 
independent stationary renewal processes each with inter-event distribution fi; we define 
also /Ug to be the probability distribution concentrated on the single-point set 5*0. For 
each n, we also regard distribution on S, assigning its total mass one to the 

set Sn- Finally, for any distribution tt on Z^, define the distribution jl^ on S hj p,^ = 
Snez+ ^iP')P'n- Thus p,^ assigns probability 7r(n) to the event that there are n individuals 
in the system, and, conditional on this event, assigns the distribution /i„ to their residual 
workloads. 

Theorem 1. Suppose that the distribution tt on Z+ is the solution of the balance equations 




'!r{n + l)/?(n -|- 1) = 7r(n)a(n), 



n > 



(1) 
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and that 

^^7r(n)a(n) < oo. (2) 

n>0 

Then the distribution p,^ on S is stationary for the process {X(t))t>o, f^nd in particular 
the distribution tt is stationary for the associated number of individuals in the system. 

Remark 1. The condition ([2|) ensures that, under stationarity, individuals arrive at the 
system at a finite rate. 

Proof of Theorem [II In order to exclude pathological behaviour in the argument below, 
we make the one additional assumption that the distribution /u has no atom of probability 
at zero. This is without loss of generality: in the case that /U does have such an atom, the 
evolution of the system may clearly be equivalently described by redefining a, (3 and // so 
as to remove it, and the result of the theorem is easily obtained via this reparametrisation. 

Analogously to the definition of for each n > 1, define the probability distribution jl^ 
on Sn to be the product of n — 1 copies of the distribution fl and a single copy of the 
distribution /i, yet again with the above identification of points in M" . (More formally, 
= A G B{Sn), where /i^") = //"~^ x fi and 9 is again the projection 

from M" into the quotient space Sn-) We again regard distribution on S, assigning 

mass one to the set Sn- 

Consider now the modified process {X{t))t>o on S describing the system in which the 
workload distribution is again /x and in which, when there are n > 1 individuals in the 
system, individual workloads are again reduced at rate /3(n)/n; however, for the modified 
system, (a) an individual departing on completion of its workload is immediately replaced 
by another bringing an independent workload with distribution /x, (b) external arrivals 
to the system are not accepted. Thus, for the modified system, the number of individu- 
als remains constant, and conditional on this being n, the system behaves as a set of n 
independent renewal processes, each of which has stationary residual workload distribu- 
tion jl. Hence, for any distribution vr' on the distribution p,^, on S is stationary for 
the process {X{t))t>o- 

Let {Pt)t>o and {Pt)t>o be the semigroups of transition kernels associated respectively with 
the processes {X{t))t>o and {X{t))t>o- For any a > 0, let Va be the class of functions / 
on S taking values in [0, 1] and satisfying the continuity condition 

\{Ptf{x) - f{x))\ <at for all X G 5 and t > 0, (3) 

where Ptf{x) = jg Pt{x,dy)f{y)- For any such / and for any distribution v on S, define 
also uf = Jg f{x)v{dx) and, for any t > 0, define vPtf = ^{Ptf) (so that fPtf is the 
expectation of f{X{t)) when {X{t))t>o is given initial distribution v); similarly define 
i^Ptf- 

Now compare the behaviour of the processes {X{t))t>o and {X(t))t>o, each started with 
the distribution p,^; so as to simplify the description below we couple these two processes 
so that they agree until the time of the first arrival or workload completion. We then have 
(see the further explanation below) that, with this common initial distribution, for any 
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o > 0, / G Va, and h > 0, 



KPhf - KPhf = Kf{X{h)) - f{X{h))) 

= hY^ 7r(n) [a(n)(/i„+i/ - fij) + (i{n){jl^_^f - fij)] + o{h) (4) 

n>0 

= ^ X] [^('^)"("') - ^("' + + 1)] (An+l/ - /^n/) + o{h) 

n>0 

= o{h) (5) 

as — > (recall in that /3(0) = 0); further the above convergence as /i — > is uniform 
over / € Pa in the sense that ^ may be written as 

sup \li^Phf - KPhf\ = o{h) as h^O (6) 

(again see below). To show @ note first that, from the above coupling and for any h > 0, 
we have f{X{h)) = f(X{h)) except where there is either at least one external arrival or at 
least one workload completion in [0, h]. It follows from the definition of fl^ that, conditional 
on the number of individuals initially being n the probability of an external arrival in [0, h] 
is a{n)h + o(h) as /i — > 0, and that an arriving individual finding the distribution of the 
system to be /2„ changes this to fj-n+i the case of the process {X{t))t>o and leaves 
it unchanged in the case of the process (X{t))t>o- Similarly, again conditional on the 
number of individuals initially being n (and recalling that m{fi) = 1), the probability of 
a workload completion in a time interval [0, h] is f3{n)h + o{h) as /i ^ 0, and that under 
the distribution conditional on such a completion taking place, the residual workload 
distribution becomes in the case of the process (X(t))t>o and fi^ in the in the case of 
the process {X{t))t>Q- Further it follows from the conditions (P) and ([2]) that, under the 
initial distribution /i^, the probability of two or more arrivals or workload completions in 
[0, h] is o{h) as /i ^ 0. That the relation ^ now holds as /i ^ with the uniformity over 
/ S T)a required for ([6]) follows easily from these results and from the definition of T>a- To 
see this note that, since / G T>a implies that / takes values in [0, 1], the contribution to 
the error term in (jH) resulting from the neglect of the possibility of two or more arrivals 
or workload completions in [0, h] is uniformly o{h) as /i — > as required. Similarly the 
terms jln+if ~ finf and fl^-if ~ f^nf ™ ([3]) are obtained by treating the precise time of 
the first arrival or workload completion within [0, /i] as if it were time h] (recalling that 
An+i' ^'^^ probability measures) it follows from ^ that the consequent error in each 
of the above two terms is bounded by 2a/i, so that the further contribution to the error 
term in ([!]) is 0{h?') as h 0, again with uniformity over / G T^a- The relations ([5]), and 
hence ([6]), are now immediate from the balance equations ([T]). Since the distribution fi^ 
is stationary for the process {X{t))t>o, it now follows from ([6]) that, again for any a > 
and h > 0, 

s^v\KPhf - Kf\ = o{h) ash^O. 
feVa 

Further, it is straightforward that if / G T>a, then also Ptf G Va for any t > 0. Standard 
manipulations using the semigroup structure of (-Pj)t>o, e.g. the consideration of increas- 
ingly refined partitions of the interval [0,t], now give that, for all a > 0, / G 2?a) and 
t > 0, 

fL^Ptf = UJ. (7) 
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Finally, we show that it follows from ([7]) that /U^i-j = fl^ for all t > 0, so that p,^ is 
stationary for {X{t))t>o as required. It is sufficient to show that, for any n > 1 and any 
set A G B{Sn) whose inverse image in M" under the mapping 9 defined above is a product 
of intervals in ]R_|_, we have 

fl^PtlA = fin^A, (8) 

where 1^ is the indicator function of the set A. It follows from the piecewise deterministic 
form of the process {X{t))t>o that we may choose a sequence of functions {fk, k > 1) such 
that, for each k, (i) fk € Va for some a > and (ii) fk and Ia agree except on a set 
whose Lebesgue measure (under 9^^) in tends to zero as A; — > oo. Since p,^, and so also 
fij^Pt, are non- atomic distributions, the result ([8]) now follows by using ([7]) with f = fk 
and letting k ^ oo. □ 

Remark 2. Suppose that the equations ([I]) above are multiplied by the signed measure 
(A„+i - fin) to give 

7r(n)a(n)(/x„+i =7r(n + l)/3(n+ -/i„), n > 0. (9) 

These equations have an obvious interpretation as representing, under the distribution p,^ 
and for each n > 0, a detailed balance of flux between 5„ and Sn+i, not just with regard 
to the total probability assigned to each of these spaces, but also with regard to the 
distribution of the residual workload sizes: the intuition underlying the derivation of ^ 
above — which is also that of Pechinkin's coupling approach — is that, under p^, an arrival 
finding n individuals in the system transforms the residual workload distribution from /2„ 
to fin+i^ while a departure from the system when it contains n + 1 individuals transforms 
the residual workload distribution from what would have been if the individual had 

remained in the system with a renewed workload, to the distribution 

In the case where we do not have processor-sharing, i.e. in which it is no longer the case 
that at any time all workloads are being reduced at the same rate, it is necessary at any 
time to distinguish the individuals in the system. Thus each Sn above is replaced by M" 
and the state space S is replaced by S* = Un>o consider as an example the case 

of the single-server queue with "last-in-first-out preemptive resume" discipline, in which 
at any time all service eff'ort is devoted to the last individual to arrive at the system. If at 
any time there are n individuals in the system, we may index these by i = 1, . . . , n in the 
order of their arrival, and no individual changes index during its time in the system; as 
usual arrivals occur as a Poisson process with rate a{n), and the workload of individual n 
is now being reduced at rate while that of the remaining individuals is being reduced 
at rate 0. As previously, define the probability distribution p on to be the stationary 
residual life distribution of the renewal process with inter-event distribution /i, and, for 
each n > 0, let the distribution p^ on be now the (ordered) product of n copies of p. 
For each n > 1, let the distribution on be the (ordered) product of n — 1 copies of 
the distribution p and a single copy of the distribution /i, with the latter assigned to the 
nth coordinate of R" . Finally, for any distribution vr on define the distribution p^ on 
S hy p^ = Yl,n&+ '^i''T-)P'n ^ before. With these (re)definitions, both Theorem [1] and its 
proof remain unchanged as stated. Again the underlying reason is as given in Remark [2] 
above: under the distribution p^, and relative to the modified process considered in the 
proof of Theorem [H an arrival finding n individuals in the system transforms the residual 
workload distribution from p^ to fL^+i^ while a departure from the system when it contains 
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n + 1 individuals transforms the residual workload distribution from An+i fin- Note that 
this balance does not obtain in the case of, for example, a "first-in-first-out" discipline, 
and here, as is again well known, we do not have the above insensitivity. 

3 Multi-class networks 

Consider now a multi-class network. We concentrate on the processor-sharing case — 
adaptations to other disciplines may be made as in the single-class case. Let I = 
{!,..., A/"} denote the set of classes, and let n = {rii, z G X) where rii is the number 
of individuals in each class i. An individual entering class i acquires a workload which has 
distribution /i* with nonzero finite mean m(/i*); we again assume without loss of generality 
that m{^^) = 1; the workload of each individual in class i is reduced at a state-dependent 
rate (pi(n)/ni, where 0i(n) > if and only if rii > 0. Individuals arrive at each class i 
from outside the network as a Poisson process with state-dependent rate <poi{n); on com- 
pletion of its workload in any class i an individual moves to class j with state-dependent 
probability (f)ij (n) / (n) or leaves the network with probability <pio (n) / <pi (n) , where 



(there are no problems in allowing the possibility 4)ii{n) > 0). The workloads, arrivals 
processes and routing decisions are all independent. 

As in the single-class case, we represent the system as a Markov process {X{t))t>o by 
defining its state at any time to be the vector n introduced above together with the 
residual workloads at that time of the set of individuals in each class. For given n, these 
workloads take values in the space Sn which is the ordered product of the spaces , i G X, 
where, as previously, each is formed from M"' by identifying points within the latter 
space which may be obtained from each other under permutation of their coordinates 
(and where again the set 5*0 contains a single point). The state space 5' for the system is 
the union of all the possible Sn, and the spaces Sn and S are endowed with the obvious 
cr-algebras B{Sn) and B{S). 

Analogously to the single-class case, for each i G Z, define the probability distribTition /7* 
on M_)_ to be the stationary residual life distribution of the renewal process with inter-event 
distribution /x' (as previously the residual life should be interpreted as a residual workload). 
For each i G X and for each > 1, define as previously the distribution Jl]^. on Sjj. to be 
the product of rii copies of the distribution /l* (again with the above identification of points 
in M"') — representing the joint distribution of the residual lives in a set of independent 
stationary renewal processes each with inter-event distribution define also JIq to be the 
probability distribution concentrated on the single-point set Sq. For each n G l/l, define 
the distribution p,^ on Sn to be the (ordered) product distribution which, for each z G X, 
assigns the distribution to . . We again regard the distribution /i„ as a distribution 
on S, assigning its total mass one to the set Sn- For any positive distribution tt on 
we define the distribution Jl^ on 5 by = X^tisz^ '^i''^)P"n- 

It is notationally convenient to expand the set X to X' = {0} U X, treating as an extra 
class feeding external arrivals to, and receiving departures from, the network. (However, 
the components of the state n of the network remain indexed in the original set X.) For 
completeness we define ^oo(»^) = for all n, and also (t>ij{n) = for all i El, j E X', and 




(10) 
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n such that = (so that pop above remains vahd for such n also). For each i (^Z, let 
Cj be the A^-dimensional vector whose ith component is 1 and whose other components 
are 0, and let bq be the A^-dimensional vector all of whose components are 0. For each n 
and each i,j G Z' define the vector T/n = n — Cj + ej; define also Tjn = T-'n = n — ei 
and T^n = T^n = n + ej. 

Theorem 2. Suppose that the distribution it on satisfies the partial balance equations 
Ti{n)Y,(t>ij{n) = Y,<T'i^)(t>ji{T',n), n € Z^', i G X', (11) 

where, for n and i ^Z such that Ui = we interpret the right side of (jlip as zero (recall 
that when Ui = we have (t>ij{n) = for all j £ Z' so that p!T]) is automatically satisfied 
in this case). Suppose also that 

^ 7r(n)^0Oi(ri) < oo. (12) 

Then Ji^ is stationary for the process [X{t))t>Q, and in particular vr is stationary for the 
associated numbers of individuals in the system. Conversely, if a distribution vr on is 
stationary for the numbers of individuals in the system for all fi = {fj,^ , i £ Z) such that 
m(;U*) = 1 for all i £Z, then vr satisfies the equations (CIP- 

Proof. Suppose first that vr satisfies the equations (llip . As in the proof of Theorem [H we 
again assume without loss of generality that each distribution /i* has no atom of probability 
at zero. 

For each n and for each i such that rij > 1 , define also the residual workload distribution 
on Sn by /i^ = fi\^ x Hj^i where /i^. is defined as in the proof of Theorem [H Thus 
corresponds to each individual in each class j having independently the stationary residual 
workload distribution /Z-', except only that a single individual in the class i is given the 
workload distribution /i*. For each n, define also fi^ = fi^. 

Again as in the proof of Theorem [H define the process {X{t))t>Q on S to be that appropri- 
ate to the modified system in which there are no arrivals, departures, or transfers between 
classes; rather each individual in each class i, on completion of its workload, acquires a 
new independent workload with distribution /i*. Thus the occupancy of the system re- 
mains constant; conditional on this being n, individual workloads in any class i such that 
Ui > 1 are again reduced at rate 4)i{n)/ni and the system behaves as a set of independent 
renewal processes. Further, for any distribution vr' on Ij^, the distribution /x^, on S is 
stationary for [X{t))t>Q. 

Again let {Pt)t>Q and (Pt)t>o be the semigroups of transition kernels associated respectively 
with the processes {X{t))t>Q and (X(t))t>o, and, for any a > 0, let Va be the class of 
functions f on S taking values in [0, 1] and satisfying the earlier continuity condition ([3]). 
Comparison of the behaviour of the processes {X{t))t>o and {X{t))t>o, each started with 
the distribution p,^ and coupled as in the earlier proof until the time of the first external 
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arrival or workload completion, now gives that, for any a > 0, / G Va, and h > 0, 

KPhf - KPhf 

= B{fiX{h))-fiX{h))) 

= ^ E ( E E ^(^)</'ji(^)AT»^/ "EE ) + o(^) 

= ^ E ( E E <nu)<\>^,{Tiu)ii^J vr(n)./.i,(n)At,/ ) + o(/i) 

= ^ E E( E^(^'^)'^i«(^'^)-E^H'^^^H)A«/+o(/i) 

„gZ^ iGX' VjgJ' jgX' / 

= o{h) (14) 
as /i — > 0, with uniformity of convergence over all / G Va, so that we may write 

sup ip^Ph/ - P-TrPhfl = o{h) as (15) 

(recall again that (f)ij{n) = whenever nj = 0, so that there is no difficulty with the 
lack of a formal definition of /x^ in this case). The identity (llSp is simply the multi-class 
version of the identity ([H) in the proof of Theorem [H and is similarly obtained, albeit 
with a slightly more compact notation; in particular, conditional on the common initial 
distribution of the two processes being given by /x^, in the time interval [0, h] where h is 
small, a transition from i to j in the original system — where either i or j may be — occurs 
with probability 7r(n)(/>jj(n)/i + o{h) as h ^ 0, and in this case the distribution of the 
process {X{t))t>o becomes /i^^ while that of the process {X{t))t>o becomes /x^; that 



()13p now holds with the required uniformity of convergence follows, using also (|12p. as in 
the earlier proof; finally the results (jl4p . and so also (jlSp . follow from the partial balance 
equations (fTTj) . Since the distribution fl^ is stationary for the process {X{t))t>o, it now 
follows from (|15p that, again for any a > and h > 0, 

sup \p-^Phf - Att/I = o{h) as h^O. 

That pL^ is now stationary for (X{t))t>o follows as in the proof of Theorem [TJ 

Now suppose that a distribution vr on is stationary for the numbers of individuals 
in the system for all /x = (/x*, i € I) with m(^*) = 1 for all i ^ Z. A proof that 
TT then necessarily satisfies the partial balance equations (jlip is given by Bonald and 
Proutiere (2002, 2003). In summary, consider the case in which in every class the workload 
distribution is exponential with mean 1, and, for any fixed class i, compare this with the 
case in which, for some < A < 1, the workload distribution in class i is replaced by 
a mixture of two distributions, obtained by choosing with probability A an exponential 
distribution with mean A~^, and with probability 1 — A the distribution concentrated on 
0. Both these models may be (re) formulated as simple Markov jump processes — in the 
latter case the transition rates into and out of the class i are reduced by a factor A. Since 
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TT is stationary in both cases, comparison of the (full) balance equations for stationarity 
yields the partial balance equations ([TT]) . □ 



Example 1. Processor-sharing Whittle networks. Suppose that for some > 0, some 
strictly positive function <I> on Z^, and some stochastic matrix P = {pij, i € 2', j G Z') 
such that Poo = 0, we have, for each n, 

(t)Qj{n) = vpQj, j^l', 

4>ij{n) = -^^Pij^ ^ ^ ^' J ^ 

where we again make the convention that ^(Tin) = whenever = 0. Then it is readily 
checked that the partial balance equations (fTTl) are satisfied by 

TTin) = a<^{n)llp^\ (16) 
for any a > and positive solution p = (pj, i E T) of the equations 

Pi = 'Y PjPji + ^Voi, i e I. 

Thus in particular the stationary distribution tt given by (1160 for the number of individuals 
of each type in the system is insensitive to the (recall our assumption m(^*) = 1 for all 
i). For the case where P is irreducible and > 0, the above equations for p have a unique 
solution. Again when P is irreducible and when 1^ = (corresponding to a closed network) 
TT remains uniquely determined, up to a multiplicative constant, by ([16]). The case where 
^{n) = Yli^x^T positive constants (Aj, i G I) characterises processor-sharing Jackson 
networks. Further discussion of Whittle networks is given by Serfozo (1999) and, for 
processor-sharing networks, by Bonald and Proutiere (2002, 2003). 

Example 2. Networks with no internal transitions. Suppose that (j)ij{n) = for all i,j € X 
and for all n € so that no transitions are possible between the classes in I. The 
partial balance equations (|lip then reduce to the detailed balance equations 

7r(n)</.io(n) = T:{Tin)<i)Qi{Tin), n G Z^", > 1, i G I. (17) 

(In the case of a single class, these equations further reduce to the equations ([I]).) An ex- 
ample is given by a traditional (uncontrolled) loss network — see, for example, Kelly (1986). 
This is naturally processor-sharing. Here workloads are identified with call durations and, 
for some set A C such that n G A implies T^n G A for all n and i such that Ui > 1 
{A is typically defined by capacity constraints), we have 

<^oi(n) = UiI{Tn G A) 
(f)ioin) = aiUi, 

for some vectors {ui, i G T) and (cxj, i G X) of strictly positive parameters. The equa- 
tions (I17p are then satisfied by 

rii 

<ri) = all^, (18) 
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where = Vi/cJi for each i, and where a is naturally chosen to be a normalising constant. 
As was originally shown by Burman et al (1984), we therefore again have insensitivity of 
the occupancy distribution vr of the network. The stationary distribution of the residual 
call durations is as identified by Theorem [2j 

Other examples of processor-sharing networks with no internal transitions are given by 
those used to model connections in communications networks with simultaneous resource 
requirements and variable bandwidth requirements — see, for example, Bonald and Mas- 
soulie (2001) and de Veciana et al (2001). Here it is far from automatic that the detailed 
balance equations (fTTl) are satisfied. 
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